Four qubits can be entangled in nine different ways. 
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We consider a single copy of a pure four-partite state of qubits and investigate its behaviour under 
the action of stochastic local quantum operations assisted by classical communication (SLOCC). 
This leads to a complete classification of all different classes of pure states of four-qubits. It is 
shown that there exist nine families of states corresponding to nine different ways of entangling 
four qubits. The states in the generic family give rise to GHZ-like entanglement. The other ones 
contain essentially 2- or 3-qubit entanglement distributed among the four parties. The concept of 
concurrence and 3-tangie is generalized to the case of mixed states of 4 qubits, giving rise to a seven 
parameter family of entanglement monotones. Finally, the SLOCC operations maximizing all these 
entanglement monotones are derived, yielding the optimal single copy distillation protocol. 
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One of the open questions in the field of quantum in- 
formation theory is to understand the different ways in 
which multipartite systems can be entangled. As the 
concept of entanglement is related to the non-local prop- 
erties of a state, local quantum operations cannot affect 
the intrinsic nature of entanglement. It is therefore nat- 
ural to define equivalence classes of states generated by 
the group of reversible SLOCC operations^ In this 
letter we are concerned with SLOCC operations on one 
copy of a state, which means that we are considering 
actions under LOCC operations on one copy of a state 
without imposing that they can be achieved with unit 
certainty. Two states belonging to the same class are 
able to perform the same QIT-tasks, although with a 
different probability. 

In the case of a single copy of an entangled pure state 
of two qubits, it is well known that it can be converted 
to the singlet state by SLOCC operations ||. In the 
case of three entangled qubits, it was shown ||, ||, |j| 
that each state can be converted by SLOCC operations 
either to the GHZ-state (|000) + |111»/V5, or to the W- 
state (|001) + |010) + |100))/V3, leading to two inequiv- 
alent ways of entangling three qubits. The GHZ-state 
is generally considered as the state with the genuine 3- 
partite entanglement, while the W-state has the pecu- 
liar property of having the maximal expected amount of 
two-partite entanglement if one party is traced out [Q. 
In this letter we extend these results to the case of four 
qubits. Furthermore the widely celebrated entanglement 
measures concurrence || and 3-tangle 0, characterizing 
the amount of genuine two- and three-qubit entangle- 
ment, are generalized to the case of four qubits, giving 
rise to a 7 parameter family of entanglement monotones. 
The SLOCC filtering operations maximizing all these en- 
tanglement monotones are derived, and it is shown that 
these are the unique operations j|] (up to local unitaries) 
bringing a state into a locally stochastic form (i.e. bring- 
ing all local density operators equal to the identity) . Fol- 



lowing Gisin |9|] , we claim that these operations maximize 
the true 4-partite entanglement. 

Interestingly, we found that there exist eight families 
of pure 4-qubit states that cannot be brought into lo- 
cal stochastic form by finite SLOCC operations. These 
states do have the peculiar property that they have 
the maximal amount of 2- and/or 3-qubit entanglement 
shared between all 4 parties. In some sense their entan- 
glement is maximally robust against the loss of one or 
two qubits. 

An interesting feature about entanglement that 
emerges out of the results of this letter is the fact that a 
quantum state has only a finite susceptibility for entan- 
glement. This will be illustrated by the fact that the op- 
erations maximizing the true 4-partite entanglement are 
precisely the operations that destroy all local correlations 
(i.e. the local density operators are made stochastic) and 
that also destroy the 3-partite entanglement (i.e. the 3- 
tanglc of the states obtained by tracing out one party 
becomes equal to zero). The states having maximal 2- 
or 3-partite entanglement shared among the four parties 
on the other hand are exactly the states having zero gen- 
uine 4-partite entanglement (i.e. the 4-concurrences are 
all equal to zero). 

Before developing the mathematical formalism, it 
should be noted that the study of states of four qubits is 
particularly interesting as the current experimental state 
of the art allows to entangle four photons |l0|, O, [l^] or 
ions |ll|. Furthermore SLOCC operations can relatively 
easily be implemented on photons, and it is therefore 
of interest to implement the optimal SLOCC operations 
such as to yield a state with maximal 4-partite entangle- 
ment. 

This letter is organized as follows. First we derive 
a simple way of determining whether two pure 4-qubit 
states are connected by local unitary operations. Next 
some advanced linear algebra is used to determine the 
orbits generated by SLOCC operations. This leads to 
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nine different families of states, corresponding to nine 
essentially different ways of entangling four qubits, al- 
though only one family is generic. This analysis gives 
rise to seven independent entanglement monotones char- 
acterizing the 4-partite entanglement. Finally the opti- 
mal SLOCC operations are derived such as to maximize 
all these entanglement monotones. 

Let us now first consider the problem to determine 
whether two pure 4-qubit states are equivalent up to local 
unitary operations. Therefore the following accident in 
Lie-group theory can be exploited: 

SU(2) ® SU{2) ~ 50(4) 

Here 15*0(4) denotes the family of real orthogonal ma- 
trices with determinant equal to 1. More specifically, it 
holds that Wx,U 2 € SU{2) : T(Ux <g> U 2 )Ti G 50(4) 
where 

(l . . 1 \ 

\t . . -ij 

A pure state of four qubits is parameterized by a four 
index tensor tpiii 2 iaii with £» G {1,2}. This tensor can 
be rewritten as a 4 x 4 matrix -0 by concatenating the 
indices (£1,12) an d (£3, £4). Next we define the matrix R 
as 

r = t$t\ (2) 

It is then straightforward to show that a local unitary 
transformation = U\ ® U 2 ® U3 ® U^\tp) results in a 
transformation R' = O1RO2 with Oi,0 2 € 50(4) and 
Ox = T(Ux <g> U 2 )T\ 2 = T(U 3 ® C/ 4 ) T Tt. A normal 
form under local unitary operations can now be imposed 
as follows: make the (1, 1) entry of R real by multiplying 
the whole matrix with the appropriate phase, and use 
Ox and O2 to diagonalize the real part of R through the 
unique real singular value decomposition. This procedure 
eliminates all 13 degrees of freedom of the local unitary 
operations, and two states are therefore equivalent up to 
local unitary operations iff they have the same normal 
form. 

Next we move to the central problem of this let- 
ter, namely characterizing the local orbits generated by 
SLOCC operations of the form 

= Ax ® A 2 <8> A 3 ® AAtj}) (3) 

with {Ai} full rank and therefore invertible 2x2 matrices. 
There is no restriction in choosing {Ai} G SL(2, C), and 
then a new useful accident arises: 

SL(2,C)®SL{2,C) ~ 50(4, C). (4) 

50(4, C) denotes the non-compact group of complex 
orthogonal matrices T = I4. Again it holds that 
VA,B G 5L(2,C) : T(A t ® A 2 )Tt G 50(4, C) with T 



given in equation (|l|), and SLOCC operations therefore 
correspond to left and right multiplication of R (|J) with 
complex orthogonal matrices. The challenge is now to 
exploit the two times 12 degrees of freedom of these com- 
plex orthogonal matrices to bring A into an unique nor- 
mal form with maximal 8 real degrees of freedom left. 
This will be possible using some advanced techniques of 
linear algebra. 

We will now state a technical theorem that is a gener- 
alization of the singular value decomposition to complex 
orthogonal matrices. 

Theorem 1 Given a complex nxn matrix R, then there 
always exist complex square orthogonal matrices Ox and 
2 such that R' = OxR0 2 is a unique direct sum of blocks 
of the form: 

1. mxm blocks of the form {Xjl m + S m ) being symmetric 
Jordan blocks (see for example \i4\l 4-4-9), and Xj is a 
complex parameter ( note that the case m — 1 corresponds 
to the scalar case). 

2. mxm blocks consisting of an upper left (mi + 1) x mi 
part being the matrix obtained by taking the odd rows 
and even columns of an (2mx + 1) x (2mi + 1) symmetric 
Jordan block, and a lower right (m — mi — 1) X (m — mi) 
part being the transpose of the matrix obtained by taking 
the odd rows and even columns of a (2(m — mi) — 1) x 
(2(m — mi) — 1) symmetric Jordan block. 

Proof: Consider the 2n x 2n complex symmetric matrix 

'-(£?)■ 151 

Due to theorem 5 in ch.XI of (IB), there exists a complex 
orthogonal Q such that P = QP'Q T with P' a direct 
sum of symmetric mxm Jordan blocks Ji with eigen- 
value Ai. Next we observe that whenever [vx \ v 2 ] (vx and 
v 2 both have n rows such that [uijua] has 2n rows) is 
the eigenspace of P corresponding to a symmetric Jor- 
dan block Ji, then [vx; —v 2 ] is the eigenspace of P corre- 
sponding to a Jordan block — Ji. Due to the uniqueness 
of the Jordan canonical decomposition, these eigenspaces 
will be either linearly independent (this holds for example 
for sure if the corresponding eigenvalue is different from 
zero) , or equal to each other (which implies that the cor- 
responding eigenvalue is equal to zero). If the first case 
applies, both vx and v 2 are orthogonal matrices. 

The second degenerated case however is more difficult. 
In this case, it holds that [i>i;v 2 ] = [vx]—v 2 ]Q for some 
orthogonal Q. Let us first calculate the standard non- 
symmetric Jordan canonical form J of the symmetric 
Jordan block with eigenvalue 0: J = UUU with U uni- 
tary and symmetric. If we define [xi;x 2 ] = [vx\v 2 ]U^ and 
Q = U^QU, the following identities hold: Q T SipQ = Sip, 
QJ = —JQ and [xx\ x 2 ] T [xx; x 2 ] — Sip (the matrix Sip is 
defined as the matrix permuting all vectors [xx, x 2 • • ■ x n ] 
to [x n , x n ~x • • • ^l])- The conditions on Q imply that Q 
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is equal to the matrix Qij 
is either of the form 



±(~l) l 8i 3 ..Therefore [a;i;x 2 ] 



Xi 

X2 



X2 



d\ bi Ci 
a 2 b 2 



ai bi 

02 62 C2 



(6) 



(7) 



Due to the constraint [xi; X2\ T [xi\ x 2 ] — Sip, the row 
dimension of [xi;x2] and therefore of J has to be odd, 
as otherwise the upper rightmost entry cannot be equal 
to 1. Rctransforming to the original picture with the 
unitary U, it holds that this structure is preserved, and 
the eigenspace [vx ; V2] is of a form t^j or ( ?]) . 

As the dimension of a Ji giving rise to the degenerated 
case has to be odd, it is compulsory that there is an even 
number of degenerated cases (indeed, the non-degenerate 
cases give rise to two times a similar block and the total 
dimension is even) . More precisely, for each [vx ; 1)2] j of 
the form @, there has to exist a [v\\V2]k of the form (Q) 
(eventually of different dimension). The eigenstructure 
of such pairs of degenerate cases can then be brought 
into the form 



a\ b\ 







a\ bf ••• 
■•■ • 



a 2 °2 





a\ b\ 



by right multiplication with a permutation matrix W. 
The effect on Ji and Jk is to transform them as 



W 



J l 
J k 



w 



( K t \ 

K k 

Kf 

V K k J 



where K u represents the matrix obtained by taking the 
odd rows and even columns of the symmetric Jordan 
block J v . 

Collecting all the pieces, it is now easily verified that 
the canonical form obtained is exactly of the form stated 
in the theorem. This completes the proof. □ 

Due to the equivalence of SL(2,C) © SL(2,C) and 
50(4, C), the normal forms arising in the above lemma 
will immediately yield a natural representative state for 
each class of 4-qubit states connected by SLOCC op- 
erations. The normal form encodes the genuine non- 
local properties of the state, while the SLOCC operators 
needed to bring the state into normal form character- 
ize the local information. The following classification is 
obtained: 

Theorem 2 A pure state of 4 qubits can, up to permu- 
tations of the qubits, be transformed into one of the fol- 
lowing 9 families of states by determinant 1 SLOCC op- 
erations (Qj: 



^(10000) + 11111)) + - 



d 



Lab 



L a b 3 — 



+£±£(|oioi) + |ioio» 

a + b 



(10000) +11111)) 



2 

a — b 



-(10110) + 11001)) 



2 

+c(|0101) + 1010)) + 0110) 
a(|0000) + 1 1111>) + &(|0101) + 1010)) 
+ 10110) + 0011) 

a + b 



(|0011) + 1 1100)) 



a(|0000) + |1111>) 
6 



+- 



10110) 



2 

1 1001)) 



(]0101> + |1010)) 



La* — 



La 2 ml 



L 



Lo 79 
3ei°3g 



+-^=(|0001) + 0010) + |0111) + |1011)) 
v2 

a(|0000) + 0101) + 1 1010) + 1 1111>) 

+(i|0001) + |0110) - i|1011)) 
a(|0000) + |1111>) + (|0011) + 0101) 
|0000) + 0101) + 1000) + 1 1110) 
0000) + |1011) + |1101) + |1110) 
|0000) + |0111) 



0110)) 



The complex parameters a,b,c,d are the unique eigenval- 
ues of P (j^j with non-negative real part, and the indices 
L a p... are representative for the Jordan block structure of 
P (e.g. ia 2 o 3e i m eans that the eigenstructure of P con- 
sists of two 2x2 Jordan blocks with eigenvalues a and 
—a, and a degenerated pair of dimension respectively 3 
and 1 ). 

Proof: If theorem is applied to a 4 x 4 R, it is easily 
checked that 12 different families arise where a family is 
defined as having Jordan and degenerated Jordan blocks 
of specific dimension. Note however that the orthogo- 
nal matrices obtained by application of the theorem can 
have determinant equal to —1, while the SLOCC oper- 
ations correspond to an orthogonal matrix with deter- 
minant +1; this is however not a problem as these op- 
erations correspond to SLOCC operations followed by a 
permutation of the qubits (1 <-> 2) or (3 4). One can 
proceed by checking that permutations of qubits (2 <-> 3) 
or (1 <-> 4) transform different families into each other. It 
is indeed true that R = Ji(a) © Ji(b) © K 3(B i transforms 
into R 1 = J2(a) © J 2 (b) if qubit 2 and 3 are permuted. 



This also happens in the case Ji(a) © K, 



501 -> J 4 (a). 
Moreover it can be shown that J\ (a) © K 3^3 is equivalent 
to Ji(a) © Js(0). Therefore only 9 essentially different 
normal forms are retained. □ 

A generic pure state of 4 qubits can always be trans- 
formed to the Gabcd state. This state is peculiar in the 
sense that all local density operators, obtained by tracing 
out all parties but one, are proportional to the identity. 
As shown in this is the unique state (up to local 
unitary operations) with this property of all states con- 
nected by SLOCC operations. In the light of the results 
of Gisin 
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we 



(10011) + 1 1100)) 



and Nielsen about majorization |L 
claim that this is the state with maximal 4-partite en- 
tanglement on the complete orbit generated by SLOCC 
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operations: the more entanglement, the more local en- 
tropy. In a later section this argument will be made hard 
by showing that a whole class of entanglement monotones 
are indeed maximized for the locally stochastic state. 

It is interesting to note that the 3-tangle JtJ of the 
mixed states obtained by tracing out one party of this 
Gated state is always equal to zero. Indeed, if the right- 
unitary matrix U 



U 



1/3 1-/3 

V2(l + |/3| 2 ) 1 ~P 1 



= y-q + V~q 2 - r 

q = 8a 2 d 2 + 8b 2 c 2 -4aV -4aV -4dV -4d 2 c 2 
r — (a 2 — d 2 )(b 2 — (?) 

is applied to the 8x2 matrix 



a + d . . a—d . b + c b—c 

b—c b+c . a—d . . a+d 



being the square root of the density operator obtained 
by tracing out the first qubit, 4 3-qubit W-states are ob- 
tained. If we define the mixed 3-tangle || as the convex 
roof of the square root of the 3-tangle, this quantity is 
clearly equal to zero. Therefore the SLOCC operations 
maximizing the 4-partite entanglement result in a loss of 
all true 3-partite entanglement. This is reminiscent to 
the case of 3 qubits where the 2-qubit state obtained by 
tracing out one particle of a GHZ-state is separable. 

Let us next discuss some specific examples. A com- 
pletely separable state belongs to the family L a &c 2 with 
a = b = c = 0. If only two qubits are entangled, an EPR 
state arises belonging to the family L a2 fc 2 w i tri a = b = 0. 
A state consisting of two EPR-pairs belongs to G a b c d 
with (a = 1; 6 = c = d = 0) or a = b = c = d, depending 
on the permutation. The class L ,„,,o*^ consists of all 
3-qubit GHZ states accompanied with a separable qubit, 
while the 3-qubit W-state belongs to the family L Q2 o 3(Sl 
with a = 0. 

The 4-qubit |$4-state belongs to the generic fam- 
ily, while the 4-qubit W-state (|0001) + |0010) + |0100> + 
|1000))/2 belongs to the family L ab3 with a = 6 = 0. This 
W-state can be shown to have a mixed 3-tangle equal to 
zero, but has a concurrence of 1/2 when whatever two 
qubits are traced out. On the contrary the state Lo 7S)i 
has all concurrences equal to zero if two qubits are traced 
out. This state is completely symmetric in the permuta- 
tion of the qubits 2,3 and 4. It has the property of having 
a mixed 3-tangle equal to 1/2 if particle 2,3 or 4 is traced 
out. This can be proven by considering the 8x2 "square 
root" 

If 1 

2 V . . . 1 . . 1 

Some straightforward calculations show that the average 
square root of the 3-tangle of the vectors obtained by 
multiplying this matrix with whatever 2 x n right-unitary 
matrix is equal to 1/2. Similar arguments show that only 
three-qubit W-type entanglement (r = 0) is retained if 
the first qubit is traced out. 



The state io Be5 is somehow a hybrid of both the 4- 
qubit W-state and Lo mi - Again a mixed 3-tangle of 1/2 
is obtained if qubit 2,3 or 4 is traced out, a mixed 3-tangle 
equal to zero if qubit 1 is traced out, but now the mixed 
state obtained by tracing out qubit 1 and (3 or 4) has a 
concurrence equal to 1/2, while the other concurrences 
vanish. 

Another interesting state belongs to the family L ai 
with a = 0: \tp) = (|0001) + |0110) + |1000>)/V3. Its 
mixed 3-tangle equals 2 /3 in the case of tracing out qubit 
1 or 4 and vanishes otherwise. Moreover the concurrence 
vanishes everywhere if 2 qubits are traced out except in 
the case of tracing out qubit 2 and 3, resulting in a con- 
currence of 2/3. 

After this zoological survey, let us next move on to the 
topic of entanglement monotones. The complex eigenval- 
ues of P (|^), given by ±(a, b, c, d), are the only invariants 
under all determinant 1 SLOCC operations (note that 
an eigenvalue is associated to the degenerated Jordan 
blocks). In H it was proven that all real positive func- 
tions of the parameters of a pure state that are linearly 
homogeneous in p and remain invariant under determi- 
nant 1 SLOCC operations, are entanglement monotones 
(in the case of mixed states they are defined by the convex 
roof formalism) . Therefore all real positive homogeneous 
functions of (a 2 , b 2 , c 2 , d 2 ) are entanglement monotones, 
such as 



M a (ip) = \a a +b a + c a + d a \ 2/a . 



Taking into account one degree of freedom due to the 
phase, this gives rise to a seven-parameter family of en- 
tanglement monotones. All these entanglement mono- 
tones are maximized by the operations making the den- 
sity matrix locally stochastic [g| (meaning that the iden- 
tity is obtained when all qubits but one are traced 
out). The optimal single-copy distillation procedure for a 
generic pure state is therefore to implement the SLOCC 
operations bringing it into its normal form G aocc i- This is 
in complete accordance with the results of Nielsen on ma- 
jorization ]lq] . Note that all the other normal forms can 
only be brought into the local stochastic normal form by 
a filtering procedure whose probability of success tends 
to zero f§. 

In summary, we have identified all different families of 
pure states of 4 qubits generated by SLOCC operations. 
Only one family is generic, and all states in it can be 
made locally stochastic by SLOCC operations. The same 
SLOCC operations represent the optimal single-copy dis- 
tillation protocol. The eight other families correspond 
to states having some kind of degenerated 4-partite en- 
tanglement and are the 4-partite generalizations of the 
3-partite W-state. 
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